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Abstract 

The elementary topological Tl process in a two-dimensional foam corresponds to the "flip" of one soap film with respect 
to the geometrical constraints. From a mechanical point of view, this Tl process is an elementary relaxation process through 
which the entire structure of an out-of-equilibrium foam evolves. The dynamics of this elementary relaxation process has been 
poorly investigated and is generally neglected during simulations of foams. We study both experimentally and theoretically the 
Tl dynamics in a dry two-dimensional foam. We show that the dynamics is controlled by the surface viscoelastic properties of 
the soap films (surface shear plus dilatational viscosity, /i s + k, and Gibbs elasticity e) , and is independent of the shear viscosity 
of the bulk liquid. Moreover, our approach illustrates that the dynamics of Tl relaxation process provides a convenient tool 
for measuring the surface rheological properties: we obtained e — 32 ± 8 mN/m and /x s + K = 1.3 ± 0.7 mPa.m.s for SDS, and 
e = 65 ± 12 mN/m and fi s + k — 31 ± 12 mPa.m.s for BSA, in good agreement with values reported in the literature. 

PACS numbers: 68.03.Cd, 68.15.+e, 83.80.Iz, 82.70.Rr 
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Foam rheology impacts material processing and prod- 
ucts in many industries and so has been the subject 
of continuous scientific activity over many years Q, 0. 
An aqueous foam acts macroscopically as a viscoelastic 
medium, whose flow depends on bulk and surface rheo- 
logical properties of the phases, which, in turn, depend on 
its constitutive ingredients (surfactant, polymers, parti- 
cles), the liquid fraction, the typical bubble size, and the 
shear rate. At low liquid fractions (a dry foam) bub- 
bles have polyhedral shapes for which local mechanical 
and thermodynamical equilibria lead to Plateau's laws 
for the main geometric characteristics: e.g., three films 
meet at each junction of a two-dimensional foam with 
equal angles of 120°. Consequently, the rheology and ge- 
ometry are linked, since as the foam structure is altered, 
rearrangements occur until a configuration is obtained 
where Plateau's laws are satisfied. 

Any rearrangement in a two-dimensional foam may be 
regarded as a combination of two elementary topological 
processes referred to as Tl and T2 pj. The Tl process 
corresponds to the "flip" of one soap film, as depicted 
in Fig. ^ while the T2 process corresponds to the dis- 
appearance of cells with three sides. From a mechanical 
point of view, the Tl process corresponds to a transi- 
tion from one metastable configuration to another, af- 
ter passing through an unstable configuration where four 
films meet at one junction (actually, for a small but fi- 
nite liquid fraction, the instability arises slightly before 
the four- fold vertex is formed Q). The spontaneous evo- 
lution from one four-fold junction to two three-fold junc- 
tions, which involves creation of a new film, is driven by 



minimization of the surface area. Various experimental 
and theoretical studies on the frequency of rearrangement 
events in foams have been conducted 0, 0, El 0] , but 
little is known about the typical relaxation time associ- 
ated with such events Indeed, the dynamics of the 
relaxation processes is usually neglected in simulations 
of foams |l(]|, UJj even though the rheological behavior of 
a foam obviously depends on this relaxation time. More 
generally, to study the evolution of the foam structure, it 
is necessary to understand the dynamics of the elemen- 
tary relaxation process. 

In this paper, we investigate theoretically and experi- 
mentally the effect of the viscoelastic parameters on the 
dynamics of the Tl process. Experiments in a two- 
dimensional foam show that the relaxation time depends 
on the interfacial viscoelasticity of the films, but not on 
the shear viscosity of the bulk liquid. These results are 
corroborated by a model, which allows for an estimation 
of the Gibbs elasticity and the surface viscosity of the 
surfactants used to make the foam. 

The experimental setup is depicted in Fig. [21 a dry 
two-dimensional foam is created in a horizontal Plexiglas 
cell (1 cm high) by blowing air through a bottle con- 
taining a surfactant solution. The polyhedral bubbles 
created have a typical edge length of 1-2 cm. The liquid 
fraction in the foam, defined as the total volume of liq- 
uids in the Plexiglas cell divided by the cell volume, is 
about 1%. Two different foaming agents have been used 
in order to study the influence of the rheological prop- 
erties of the interface on the Tl dynamics: (i) Sodium 
Dodecyl Sulfate (SDS), at a concentration of 4.80 g/L, 
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FIG. 1: (color online) Schematic of the Tl transition. The ini- 
tial configuration (a) evolves continuously through metastable 
states, for which Plateau's laws are satisfied, to an unstable 
four- fold configuration (b). This unstable state spontaneously 
evolves into two three-fold junctions with creation of a new 
film (c) until a new metastable configuration (d) is reached, 
and Plateau's laws are satisfied again. Topologically, the tran- 
sition between the initial configuration (a) and the final con- 
figuration (b) corresponds to a "flip" of one soap film. 



forms "mobile" surfaces, and (ii) protein Bovine Serum 
Albumin, together with a cosurfactant Propylene Glycol 
Alginate (PGA), both at concentrations of 4.00 g/L, form 
"rigid" interfaces. The impact of the shear viscosity of 
the bulk liquid has been investigated by adding glycerol, 
0%, 60% and 72% (w/w), to the bulk solutions. 

In order to cause rearrangements in the foam, we use a 
syringe to blow the air away from one bubble. Then the 
rearrangements are viewed from above with a high-speed 
camera. The length of the soap film is measured by fol- 
lowing its two ends using particle-tracking software. Ex- 
periments where the other vertices have noticeable move- 
ments during the relaxation process, or where simultane- 
ous Tl events occur, have been disregarded. 

The relaxation process is characterized by the creation 
of a new soap film following the appearance of a four-fold 
junction. For different surfactant systems, we report the 
length of the new film, normalized by its final length, as 
a function of time in Fig. Several trials for each so- 
lution are shown and illustrate that the time evolution 
of the reduced length appears independent of the final 
film length. We define a typical time T associated with 
the relaxation process as the time for the film to reach 
90% of its final length. A comparison of the results for a 
foam made with the SDS without glycerol (shear viscos- 
ity of the liquid [i — 1.0 mPa.s) and with 60% glycerol 
(/i = 10.7 mPa.s) shows that there is no significant effect 
of the viscosity of the bulk liquid (T ~ 0.5 sec for both so- 
lutions). This response is not unreasonable since, for the 
mechanics of a free soap film, viscous effects of the bulk 
are generally negligible in comparison with the effects of 




FIG. 2: (color online) Experimental setup: a two-dimensional 
foam is created in a 1 cm high horizontal Plexiglas cell by 
blowing air through a surfactant solution. The bubbles are 
polyhedral with an edge length of 1-2 cm and the liquid frac- 
tion in the foam is about 1%. Using a syringe, air is blown 
away from one bubble, which induces rearrangements, that 
are viewed from above using a high-speed camera. 



the viscoelastic properties of the interfaces [12J, [13| . 

We note that as the foam is constrained between two 
planes, most of the liquid of each soap film is located 
in the menisci close to the solid surfaces. In this region 
there is dissipation that depends on the shear viscosity of 
the solution. Hence, from our experimental observations 
we conclude that frictional effects at the boundaries have 
a negligible influence on the Tl dynamics. This result is 
not in contradiction with the observations made on the 
rheology of 2D foam 0, 0, 0] , where a macroscopic 
stress (pressure drop) causes motion of the foam relative 
to the boundaries and the shear viscosity of the bulk 
solution is the main parameter controlling the dynamics. 

Next, we compare the SDS results with those from the 
BS A/PGA solution /j, = 7 mPa.s, both without glyc- 
erol. We observe a change in the typical relaxation time 
by about a factor of 7 (T ~ 3.7 sec for the BSA/PGA 
solution). Although, coincidentally, the viscosity of the 
BSA/PGA solution is increased by a factor 7 relative to 
the SDS solution, we can rule out this influence since we 
just demonstrated that the viscosity of the solution is 
not rate limiting. We also verified that the addition of 
glycerol to the BSA/PGA solution did not produce any 
significant change in the typical relaxation times (results 
not shown). These results allow us to conclude that the 
viscoelastic properties of the interfaces dictate the relax- 
ation time of the Tl process. 

We now provide a brief description of a theoretical 
framework for the Tl dynamics; for details of the deriva- 
tion and complete considerations of various special cases 
see |l7J • We compare the theoretical predictions with the 
experiments, which allows us to extract surface Theolog- 
ical parameters. The net result is a model for the Tl 
process and estimates for the relaxation time. 

We assume that the geometry is symmetrical with thin 
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FIG. 3: (color online) Evolution of the film length, normal- 
ized by its final length, with time. Each curve represents 
one experiment. Red curves: SDS solution without glycerol; 
green curves: SDS solution with 60% (w/w) glycerol; blue 
curves: BS A/PGA solution without glycerol. The data for 
SDS solutions with glycerol overlay those without glycerol, 
which confirms that viscosity of the bulk liquid is not signif- 
icant. The typical time of the relaxation process, defined as 
the time to reach 90% of the final length, is about 0.5 sec for 
the SDS curves and 3.7 sec for the BSA/PGA curves. 



films nearly together at an unstable four-fold configura- 
tion; each film connects to one of four fixed vertices at 
the corner of a rectangle with sides 2L X and 2L y (see Fig. 
^) . It is experimentally observed that the stretching film 
has a spatially uniform but time varying thickness h(t) 
everywhere except near the Plateau borders located at 
x = ±XB{t) at the top and bottom boundaries; these dy- 
namics are common for the fluid dynamics of thin films 
flif . We focus on the dynamics of the stretching central 
film of length 2xb (t) , which is driven by the monotonic 
decrease of the angle a(t): 



cosa(i) = 



L x - x B {t) 



yjLl + (L x -x B (t)Y 



(1) 



The surface tension 7 changes in time as the surfac- 
tant surface density T is reduced by stretching. Also, we 

introduce the length L(t) = ' L? y + (L x — x B (t)) 2 of the 
adjacent films which shorten with time. Since the flat 
film can support no pressure gradient the axial velocity 
U(x, t) in the main body of the film is a linear function 
of position, U oc x This uniform extensional motion 
means that away from other films we have ^(t) and 
which do not change with position x. 

The dynamics follow from a force balance on the 
stretching film and a surfactant mass balance. Neglect- 



ing inertia terms and dissipative terms associated with 
the viscosity of the bulk liquid 0], Newton's second law 
applied to the film reduces to a balance between surface 
tension contributions and surface dissipative terms (fi s 
and k denote the shear and dilatational viscosities and 
7e 9 denotes the equilibrium value of the surface tension 
of the soap solution): 

dlf 

2j eq cos a(t) - 7 (i) - (fi s + K ) — = 0. (2) 

During the expansion of the new film, adjacent films act 
as surfactant reservoirs, so that j(t) and T(t) are as- 
sumed to be close to their equilibrium values, and so 
are related by the Langmuir equation of state: j(t) = 
7 eg — eln(r(i)/F e9 ), where e is the Gibbs elasticity and 
T eq is the equilibrium surface density. This equilibrium 
value T eq is assumed to be present in the adjacent films. 

Next, we turn to a mass balance on the surfactant. 
This requires accounting for stretching of the interface 
as well as addition of surfactant to the new surface cre- 
ated as the adjacent films, of length L(t), are shortened. 
We assume that as surfactant is added to the end of the 
stretching film the surface density is rapidly adjusted to 
a spatially uniform state by strong Marangoni forces. In 
addition, we neglect diffusion/adsorption processes from 
the bulk liquid, restricting our study at short times (these 
slower processes only affect shape adjustments as the fi- 
nal equilibrium is approached). So, the change in the to- 
tal number of surfactant molecules during a time interval 
dt is d(Txb) = —T eq dL, where dL — —dx B .cosa(t); here 
we have only accounted for one interface of each adjacent 
film as feeding the stretched film since surfactant from 
the other interface must desorb from the surface, tran- 
sit through the bulk, then adsorb to the stretched inter- 
face, which is a much longer process. Integration of this 
equation allows us to express T(t) as a function of x B (t): 

-L(t 



T(t) 

(xo - 



= r, 



where L c = \ IA + (L x - x ) 



e 1 x B (t) ' -"C — y ~y 

^b(O)). In addition, with the above approxima- 
tions, the density along the film must also satisfy the local 
conservation law ^ +T^r = 0. Comparison of these two 

evolution equations leads to ^ = 2s. M _ cos a(t)j . 
Note that since U cx x, the surface velocity at the junc- 
tion is U(x B (t),t) = (l - §f^cosa(t)\ which is 

smaller than the velocity of the junction itself (= x B )- 
This velocity difference is a consequence of the slip of 
the surface (and surfactant) coming from the adjacent 
film. Finally, Eq. [21 is rewritten as an evolution equation 
for x B (t): 



2-f eq ( cos a(t) - i 



In 



(Ms + K ) 



x B 
x B 



1 - 



L c - L{t) 
x B (t) 

x B (t) 



L(t) 



cosa(i) = 0. (3) 
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We now compare this theoretical description with the 
experiments using measurements of xb ft) ■ From Eq. [31 

a plot of Y(t) = If- (l — cos aft)) / (cos aft) - \) 

versus Aft) = In ( L ^2^j^ J I (cos aft) — |) should yield 

-X. 



XB(t) 

a straight line: Y = 2 



Typical experimental curves Y(t ) vs Aft) are shown 
in Fig. H for SDS and BSA/PGA |2l|. At short times 
(typically, for times below 0.1 sec for SDS and 0.8 sec 
for BSA), their evolution is nearly linear, in excellent 
agreement with the theory. A linear fit of the experi- 
mental data plotted in this way allows determination of 
fi s + k and e, using well-established equilibrium surface 
tension values of SDS and BSA (38 mN/m and 55 
mN/m |23| . respectively). The mean values we obtained 
are e = 32 ± 8 mN/m and fi s + K = 1.3 ± 0.7 mPa.m.s for 
SDS, and e = 65±12 mN/m and fi s +n = 31±12 mPa.m.s 
for BSA/PGA, which are in good agreement with values 
reported in the literature (e.g. [3 EH)- Note, however, 
that our value of the surface viscosity of SDS correspond 
to the highest value reported in literature [24|, which, in 
both studies, may be a consequence of the rapid surface 
stretching giving rise to some nonlinear or inertial effects. 
In addition, there is a large difference between values of 
the shear and dilatational surface viscosities, and it is of- 
ten unclear in published works which surface viscosity is 
actually measured. 

This study of the Tl dynamics in a two-dimensional 
foam illustrates that the relaxation time T associated 
with the process is a function of two parameters, Ps+K 

and < J,a J K . From dimensional analysis of Eq. |3 
A ' 3+K /( — ), where / is an increasing function of the di- 

7eg Teg 

mensionless parameter — . This theoretical description, 
which has been corroborated with our experimental data, 
might be useful for simulations of aging or rheological 
properties of foams. Finally, we expect that a sheared 
foam has a different rheological response when the shear 
rate is significantly different than this typical relaxation 
time. 
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